Three-body amplification of photon heat tunneling 
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Resonant tunneling of surface polaritons across a subwavelength vacuum gap between two polar 
or metallic bodies at different temperatures leads to an almost monochromatic heat transfer which 
can exceed by several orders of magnitude the far-field upper limit predicted by Planck's blackbody 
theory. However, despite its strong magnitude, this transfer is very far from the maximum theoretical 
limit predicted in near-field. Here we propose an amplifier for the photon heat tunneling based on a 
passive relay system intercalated between the two bodies, which is able to partially compensate the 
intrinsic exponential damping of energy transmission probability thanks to three-body interaction 
mechanisms. As an immediate corollary, we show that the exalted transfer observed in near-field 
between two media can be exported at larger separation distances using such a relay. Photon heat 
tunneling assisted by three-body interactions enables novel applications for thermal management at 
nanoscale, near-field energy conversion and infrared spectroscopy. 



I. INTRODUCTION 

Two bodies held at different temperatures and sep- 
arated by a vacuum gap exchange in permanence heat 
throughout the thermal electromagnetic field they ra- 
diate. At long separation distance (i.e. in far field) 
this exchange of energy results exclusively from prop- 
agative photons emitted by these sources. The black- 
body limit given by the famous Stefan-Boltzmann's law 
sets the maximum heat flux these media can exchange. 
However, at subwavelength distances (i.e. in near-field 
regime) the situation radically changes. Indeed, at this 
scale, evanescent photons which remain confined near the 
surface of materials are the main contributors to transfer 
and they participate via tunneling through the vacuum 
gap [THE]. A significant heat-flux increase results from 
this transport. In presence of resonant surface modes 
such as surface plasmons or surface polaritons, collective 
electron or partial charge oscillations coupled to light 
waves at the surface, the radiative-heat exchanges can 
even drastically surpass by several orders of magnitude 
the prediction of Planck's blackbody theory. This discov- 
ery opened new possibilities for the development of inno- 
vative technologies for nanoscale thermal management, 
such as near-field energy conversion (thermophotovoltaic 
conversion devices [3110]), heating assisted data storage 
(plasmon assisted nanophotolitography [11]) or IR sens- 
ing and spectroscopy [HI [T3] . 

Recent theoretical developments [21 [TS] have estab- 
lished that the nanoscale heat transfer between two me- 
dia separated by a distance d can be revisited using the 
same concepts as in mesoscopic physics of charge trans- 
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port. The first step is the expression of the distance- 
dependent heat flux ip(d) under the form of a spec- 
tral decomposition <p(d) — J Q °° ^(f>(uj,d). Hence, the 
monochromatic near-field heat flux transferred between 
two media at temperatures T\ and T 2 is described in 
a Landauer-like formalism [TIHT7] as the sum over all 
modes 

p * ' 

of the energy hui of each mode (u;, k,f>) (identified by 
the frequency u>, the transverse wavevector k = (k x ,k y ) 
and the polarization p), times a transmission probabil- 
ity T p (uj, k, d) through the separation gap (assuming val- 
ues between and 1), with nij(ui) = n(uj,Ti) — n(to,Tj), 
n(uj,T) = (e hM f kBT — l)' 1 being the distribution function 
inside the reservoir of modes at temperature T. Analy- 
sis of arbitrary systems has revealed a rather small ef- 
ficiency of transfer coefficient for each exchange channel 
even in the presence of resonant surface modes such as 
surface polaritons. In this situation the total number 
of channels (i.e. number of modes) contributing to the 
heat transfer per unit of surface is set by the position 
of the cutoff for the transverse component of wavevec- 
tor k c ~ log[2/ v / Im(e 1 )lm(e 2 )]/d [HJOI], where Im ( e i) 
denotes the imaginary part of the dielectric permittivity 
of the material i (i — 1,2). This allows to estimate two 
fundamental limits [14] max = kg(Tf — T|)fc 2 /6fi. and 
/i max = <7o^c/ 7r (.9o — Tr 2 kgT/3h is the quantum of ther- 
mal conductance at temperature T) respectively for the 
monochromatic flux and total heat transfer coefficient be- 
tween two media. However, these limits are intrinsic to 
two-body systems but are not universals in general. In 
the present work we describe heat transfer in three-body 
systems and we highlight the concept of three-body am- 
plification of heat flux exchanged at nanoscale between 
two media. In addition we propose a device based on an 
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intermediate passive relay which is able to increase the 
number of coupled modes. This investigation belongs to 
the vast category of few-body problems whose richness 
has been largely explored in atomic physics, quantum 
chemistry and celestial mechanics. Recently, interesting 
effects have been theoretically discussed in the context of 
heat transfer [11120]. 



II. RESULTS 

A. Physical system and analytical results 

We consider a system consisting of two parallel slabs 
identified for convenience by the indexes i = 1 and i = 3, 
as shown in figure [ija). Each slab has a finite thickness, 
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FIG. 1: Geometry of the system in the (a) two- and (b) three- 
slab configurations. The distance d between adjacent slabs 
(1-2 and 2-3) in the three-body configuration is chosen to be 
equal to the distance between slabs in the two-body case. 

whereas its transverse extension is much larger than the 
distances between the slabs, so that the planar surfaces 
can be considered infinite with respect to the xy plane. 
We compare this system to a configuration in which a 
third slab, labeled with i = 2 and having thickness 5, 
is placed between slabs 1 and 3 (see figure [ljb)). The 
system is placed in both cases in a stationary thermody- 
namical configuration, in which each body is held at a 
temperature Ti (i = 1,2,3). Moreover, in order to take 
into account the radiation impinging on the system from 
outside, we assume that the two slabs are embedded in 
an environment described as a source of isotropic ther- 
mal radiation at a constant temperature T e , in general 
different from the temperatures Ti. As far as the material 
properties of the three slabs are concerned, we describe 
them in terms of three complex dielectric permittivities 
£i(uj), meaning that our bodies are characterized here by 



a non-magnetic local electromagnetic response. 

Herein we are interested in comparing the heat flux 
on body 3 in two- and three-body configurations. The 
heat-transfer problem in the case of a couple of arbitrary 
bodies in a thermal environment has been recently solved 
|21j . In this case, the heat flux on a body can be ex- 
pressed as a sum of an evanescent and a propagative con- 
tribution: the former depends only on the modes of the 
electromagnetic field for which the transverse wavevector 
k satisfies ck > uj (i.e. outside the light cone), whilst the 
latter is defined in the region ck < lu. While the propaga- 
tive term depends on all the three temperatures involved 
in the system, the evanescent contribution is only con- 
nected to the disequilibrium between the temperatures 
Ti and T 3 of the two bodies. Then, for T x ^ T 3 , at 
distances between them smaller than the thermal wave- 
length (of the order of some microns at ambient tem- 
perature) the evanescent contribution largely dominates 
on the propagative term. In this near-field regime the 
monochromatic heat flux on body 3 at frequency oj can 
be written under the form of a Landauer expansion 



^ni 3 (w)V / j— ^T 2s ,p(a;,k,d) (2) 
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These quantities depend also on the thicknesses of slabs 
1 and 3 through the reflection coefficients pi of slab i [22] . 

The scattering procedure developed to investigate heat 
and momentum transfer between two bodies described in 
detail in [22] has been generalized to the case of three 
bodies in a thermal environment [23]. In this case, as 
expected on physical grounds, the evanescent contribu- 
tion is a function of the three temperatures Ti, T 2 and 
X3, while the environmental temperature T e enters again 
only in the propagative term. In order to reduce the 
number of free parameters for the analysis of the am- 
plification mechanism, we will focus our attention of the 
symmetric case in which the distances between adjacent 
slabs in the three-body configuration (1-2 and 2-3) are 
both equal to d, i.e. the distance between slabs 1 and 3 
in the two-body configuration (see figure IJ. This choice 
makes d the only relevant distance in both scenarios. 
Moreover, it makes the minimal distance between any 
couple of adjacent bodies the same for the two configu- 
rations and the optical distance between slabs 1 and 3 
double in the three-slab case with respect to two slabs. 
For the three-slab system the near-field expression of the 
monochromatic heat flux on slab 3 takes the form 123 



b 3s (uj,d,S) = 



d 2 k 

(23), 



n 12 (u)T^(u},^d,5) 



+ «23(w)73s jP ; (w,k,d,(5) 
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with 
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Equation Q shows that the evanescent term splits into 
two contributions, associated to the couples of tempera- 
tures (T\,T2) and (T2,Ts), each one containing a differ- 
ent Landauer transmission probability which is weighted 
using the corresponding difference (w) of thermal pop- 
ulations. In this case, the dependence of the transmission 
amplitudes on the thickness S of the intermediate slab has 
been explicitly stated. The amplitudes ^ and (6j) are 
functions of the reflection and transmission amplitudes 
associated to each slab and to the couple (1,2) treated as 
a unique body, taking the form 



Pi2 P (k,w) = p 2p (k, u) 



r 2p( k »^)Pip( k ^) e 
|1 - pip(k,w)p2 P (k 
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Let us now consider the choice of temperatures. To 
avoid the introduction of a supplementary thermal source 
and thus keep the comparison between the two configu- 
rations meaningful we choose for the temperature T2 the 
value such that the heat flux on body 2 is zero. In the par- 
ticular case of a quasi-monochromatic flux at frequency 
ujq (which, as we will see, is quite accurate in near-field 
regime) it can be shown that T2 is solution of the equa- 
tion 2n(wo,T2) = n(uJo,Ti) + n(uiQ,T 3 ). We are finally 
left with the discussion of the material properties of the 
three slabs. Regarding media 1 and 3, we have chosen 
two SiC slabs of thickness 5 pm. As for the intermediate 
slab, it is a metallic-like medium described by a simple 
Drude model s(uj) = 1 — u>p/uj(u! + ijp), defined by a 
plasma frequency up and a relaxation rate jp. As it is 
well known [24], this model predicts the existence of a 
surface mode, a plasmon, at a frequency close to cjp/%/2- 
Hence, a natural choice for wp to maximize heat transfer 
is ujp — \/2w S pp, where ui spp ~ 1.787 ■ 10 14 rads~ 1 is the 
surface plasmon-polariton frequency supported by the 
SiC slabs, matching the SiC and Drude surface-polariton 
frequencies. 



B. Heat-flux amplification 

We are now ready to calculate the ratio between the 
total heat flux on slab 3 in the three-slab case with re- 
spect to the one in presence of two slabs only. Figure 
[2] shows the value of this ratio in a range of distances 
d £ [50, 800] nm and thicknesses S € [0, 1] pm of the cen- 
tral slab. This figure clearly shows the existence of a 
region of parameters in which the three-body heat flux is 
larger than its two-body counterpart. As far as the value 
of the amplification is concerned, in our physical config- 
uration it goes up to 60% for short distances (around 




FIG. 2: Total heat-flux amplification ips s (d, 8)/tp2s(d) as a 
function of distance and slab thickness. The black dashed 
line corresponds to v?3 S (d, 8)/ip2s(d) = 1. The dielectric per- 
mittivity of SiC is described using a Drude-Lorentz model 
[25] , and for slab 2 we have chosen jp — lCP 3 ajp. The inset 
shows the dependence of amplification rates on the relaxation 
rate 7p for d = 200nm (black solid line), d = 400 nm (green 
dashed line) and d — 600 nm (blue dot-dashed line). For each 
value of d the associated optimal value of 8 has been used. 



150 nm). As shown in the inset of figure [2j the amplifica- 
tion effect is robust even with respect to the relaxation 
rate jp over a range going up to values typical of met- 
als. Moreover, the amplification exists up to distances 
around 650 nm and, for each given value of the distance 
d, is limited to a finite interval of 6. More in detail, for 
a given value of d, we see that the amplification assumes 
a maximum as a function of the thickness of slab 2 for 
5~d. 



III. DISCUSSION 

In order to get an insight on the highlighted ampli- 
fication mechanism, we first analyze the flux spectra 
</>2s(w, d) and 4>3 S (uj,d,8) for different values of the dis- 
tance d and for the optimal value of the thickness 6 as- 
sociated to each value of d. The results plotted in fig- 
ure for d = 200 nm and d = 700 nm show that when 
the distance between the two slabs approaches zero, the 
flux tends to become monochromatic at the frequency of 
the surface polariton. As far as the three-slab term is 
concerned, we see that in both cases the spectral com- 
ponent of the flux at w spp is enhanced with respect to 
the two-body case. Nevertheless, the presence of slab 2 
does not amplify the contribution at smaller frequency, 
but on the contrary provides a quasi-monochromatic en- 
hancement around the surface polariton frequency o; spp . 
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This explains why, when increasing the distance d, the 
amplification is less important and eventually the ratio 
between the fluxes goes below one. 
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FIG. 3: Monochromatic fluxes c/>2 S (^, d) (red dashed line) and 
c/>3 a (aj, d, 8) (black solid line) in arbitrary units for two differ- 
ent values of the distance d. For each distance the associated 
optimal value of 8 has been used. 

This behaviour can be understood by looking at the 
coupling efficiencies of modes defined in equations ([3]), 
|5]) and (J6|. It is important to remark that in the two- 
and three-body configurations the transmission probabil- 
ities are multiplied by different weighting functions n^. 
In the two-body case, this weight is only given by the 
difference 7113, while 7742 and n-ij, appear in the case of 
three slabs. However, due to the quasi-monochromaticity 
of the flux, the temperature T2 is such that ni2(w S pp) = 
^23(w spp ) = "-i3(w spp )/2. Hence, we compare the trans- 
mission probability T 2s to T 3s = (T 3 [ 12) +T 3 ( s 23) )/2. These 
quantities are plotted in figure [4] in the plane (w, k) for a 
given couple (d, 8). As a first observation, we see that the 
resonant peak in the spectral flux plotted in figure [3] is 
due to the presence of resonant surface modes which exist 
only in TM polarization. We then note that, while in the 
two-body case we have one symmetric and one antisym- 
metric surface mode [21], we observe the appearance of 
two supplementary modes in presence of an intermediate 
slab which results from doubling the number of cavities 
in the system. Moreover, as shown in the right side of 
figure |4j in the three-slab case these modes remain effi- 
cient up to larger values of k. These features result in 
the enhancement and broadening of the resonance peak 
of the spectral flux. For both polarizations, the presence 
of supplementary (frustrated) modes at low frequencies 
explains the secondary peak in the spectral flux in figure 
[3j The number and efficiency of those modes is reduced in 
the three-body configuration: this explains why the am- 
plification mechanism is centered around the resonance 
frequency of both surface modes (figure [3]). 

To go forward in the physical analysis of the exaltation 
mechanism we now focus our attention on the transmis- 



sion probabilities 72s and 73 S (figure [5]) at the resonance 
frequency w spp as a function of the wavevector k and 
thickness S for a separation distance d = 200 nm. We re- 
strict ourselves to TM polarization, which gives the main 
contribution to the heat transfer. In figures]^ a) and[5jb) 
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FIG. 4: Transmission probabilities 72s (upper side) and TL 
(lower side) for both polarizations as a function of lu and k at 
separation distance d — 700 nm and for the associated optimal 
value of 8. 

we show the two- and three-slab transmission probabili- 
ties 72s and 75s respectively. Let us start by discussing in 
figure^b) the two limiting regimes with respect to 8. For 
8 = the system reduces to two SiC slabs separated by a 
distance 2d. By comparing this case to the two-slab con- 
figuration at distance d in figure we observe a shift 
of the peak and of the cutoff wavevector k c correspond- 
ing to the difference of distance from d to 2d. Otherwise, 
when the thickness 8 is larger than d, Tz s becomes inde- 
pendent on 8, corresponding to the fact that the slab 2 is 
seen as semi-infinite. In this limit, figure^b) shows that 
the replacement of the SiC slab 1 with a semi-infinite slab 
described by the Drude model produces a shift of k c to- 
ward larger values. This is consistent with the fact that 
k c is an increasing function of 1/ -\/lni(ei)Im(e 2 ) and that 
at the plasmon frequency uj spp the imaginary part of the 
Drude dielectric permittivity is smaller than the one of 
the Drude-Lorentz model describing the SiC. Despite this 
shift, the integral over k (i.e. the monochromatic flux) is 
still smaller than in the case of two slabs at distance d. 
This is a direct consequence of the fact that T2 < Ti. 

Hence, we can say that the cases 8 = and 8 3> d are 
as a matter of fact two-body configurations (SiC - SiC and 
Drude -SiC respectively). The intermediate region, cor- 
responding to values of 8 around d, results from a purely 
three-body effects in which the results is intimately con- 
nected to the presence of both bodies 1 and 2. A deeper 
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understanding of this transition comes from figures [5jc) 
and[5jd), representing separately the two contributions 



to the three-slab configuration 0.57^ 



(23) 
3s 



and 0.5T 3 ( S 12) 



spectively. The former has exactly the structure of a two- 
body transmission amplitude, as evident from the com- 
parison between equations ([3| and ([6]). It corresponds to 
the exchange between the couple (1,2) treated as a unique 
body at temperature T 2 and the body 3. The plot of this 
coefficient in figure [5^c) shows that this term can be asso- 
ciated to the transition, as a function of 5, from body 2 to 
body 1. As a matter of fact, the peak of the transmission 
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FIG. 5: (a) Transmission probability Tis = 0.5(7^' s 12 ' ) + 
7^' s 23 ') of a three-body SiC-Drude-SiC system at the sur- 
face plasmon-polariton frequency Lo S p P and at a separation 
distance d — 200 nm as a function of dimensionless wavevec- 
tor ck/uj sp p and thickness 5; (b) 72 S ; (c) 0.57Jjg 23 '; (d) 0.5T^ 2 ^ ; 
(e) Particular case of optimal thickness S = 265 nm. 



temperature-dependent term ni2(w) and in the fact that 
the product of imaginary part of reflection coefficients p\ 
and p 3 is weighted over a coefficient dependent on body 2 
and in particular proportional to the square modulus of 
its transmission coefficient. This explains why this term 
starts contributing below a given value of 5 (of the order 
of 300 nm in figure), at which the slab 2 is no longer seen 
as semi-infinite. It is in proximity of this value that the 
system exploits the shift of the cutoff wavevector due to 
the replacement of the SiC with the Drude material and 
at the same time the contribution for smaller wavevec- 
tors guaranteed by the presence of slab 1. This is even 
more evident in figure [5je) . In this case we represent the 
two-slab configuration at distance d, and the three-slab 
cases for 5 = 0, S ^> d and i5 = 265 nm, representing in 
this case the optimal value. The curve corresponding to 
this value of 5 shares the behaviour for k — ¥ oo with the 
case S ^> d, but differs from this curve at smaller values 
of k. 

In conclusion, we have proposed and characterized a 
new passive amplifier based on a three-body assisted 
tunneling mechanism. By remarkably enhancing near- 
field exchanges with respect to a two-body system, it 
allows at the same time to increase the magnitude of 
heat transfer and to double the separation distance, with- 
out any additional source of energy. Due to its quasi- 
monochromaticity and filtering effect, this mechanism 
could also be exploited to improve considerably the effi- 
ciency of near-field energy conversion devices by increas- 
ing the photocurrent generation in thcrmophotovoltaic 
cells, as well as reducing the heating of the cell. This work 
also paves the way to the study of near-field heat trans- 
port in complex plasmonic systems mediated by many- 
body interactions at mesoscopic scale. 



amplitude is for any thickness approximately 0.5, but the 
range of wavevectors (i.e. the number of modes) con- 
tributing to the effect moves from the one associated to 
a Drude material at distance d to the one corresponding 
to a SiC slab at distance 2d. On the contrary, the con- 
tribution 7^g 12 ' can be thought as an exchange between 
bodies 1 and 3 mediated by the presence of body 2. The 
presence of the intermediate slab is manifest both in the 
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